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This paper presents an experimental and theoretical study of the distri- 
bution of carbon atoms in the octahedral interstitial sites of the face-centered 
cubic {fee) phase of the iron-carbon system. The experimental part of the 
work consists of Mossbauer measurements in Fe-C alloys with up to about 12 
atomic percent C, which are interpreted in terms of two alternative models 
for the distribution of C atoms in the interstitial sites. The theoretical part 
combines an analysis of the chemical potential of C based on the quasichem- 
ical approximation to the statistical mechanics of interstitial solutions, with 
three-dimensional Monte Carlo simulations. The latter were performed by 
assuming a gas like mixture of C atoms and vacancies (Va) in the octahedral 
interstitial sites. The number of C-C, C-Va and Va-Va pairs calculated using 
Monte Carlo simulations are compared with those given by the quasichemical 
model. Furthermore, the relative fraction of the various Fe environments were 
calculated and compared with those extracted from the Mossbauer spectra. 
The simulations reproduce remarkably well the relative fractions obtained 
assuming the FegCi-y model for Mossbauer spectra, which includes some 
blocking of the nearest neighbour interstitial sites by a C atom. With the 



new experimental and theoretical information obtained in the present study, 
a critical discussion is reported of the extent to which such blocking effect is 
accounted for in current thermodynamic models of the Fe-C fee phase. 
PACS Codes: 2.70.Uu, 76. 

I. INTRODUCTION 

The physical properties of the austenite solid solution phase have been studied extensively 
over the years in connection with, e.g., the assessment and understanding of the phase 
diagram [1-6], the diffusion controlled [7] and martensitic phase transitions in Fe-C alloys 
[8]. In austenite the iron atoms are arranged in a close-packed face-centred cubic {fee) lattice, 
and the C atoms occupy a limited number of the octahedral interstices which are located at 
the centres and at the mid-points of the edges of the unit cubes, these two positions being 
crystallographically equivalent (Fig.l) [9]. Accordingly, various models of austenite have 
been proposed which are based on assuming two sublattices, one for the Fe atoms and the 
second one for the mixture of carbon atoms (C) and vacant octahedral interstices (Va). The 
general theme of the present paper is the distribution of the C atoms in the interstitial sites, 
as revealed by three complementary sources of information, viz., thermodynamic properties, 
Mossbauer experiments and Monte Carlo simulations. 

In the ideal solution model (ISM) for thermodynamic properties of austenite it is assumed 
that C and Va distribute themselves at random in the octahedral interstitial sites, the amount 
of which is equal to the number of Fe atoms. On these basis, the thermodynamic activity 
(oc) of C in an ideal mixture of carbon atoms with Npg iron atoms is shown to be 
proportional to the ratio yc/(l — He), where i/c — Nc/N^e, and 1 — i/c represent the fraction 
of occupied and of empty interstitial sites, respectively [10-13]. Since the experimental 
Oc in austenite deviates positively from the ISM, many approaches have been proposed 
to account for the non-ideal behaviour [14-22]. The reader is referred to ref. [23] for a 
recent review of the work of most relevance for the present study. In the strict version 
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of the approach known as the hard-blocking excluded-sites model (HBESM) it is assumed 
that the presence of a solute atom blocks the occupancy of a certain number (b) of the 
nearest neighbor interstitial sites (NNIS), so that a site is either blocked or is available 
for the mixing of C atoms and Va [9, 15 , 23]. Further, if the mixing in the non- blocked 
sites occurs at random, the Uc in austenite becomes proportional to yc/[l ~ (1 + b)yc] [23]. 
Frequently, b has been treated as an adjustable parameter, identified with the value to be 
inserted in the expression for Oc in order to reproduce the experimental data [9,19,22,23]. 
Alternatively, some theoretical studies have been reported which suggest that b should in 
fact be treated as composition dependent [18,19]. A different approach will be explored 
in the present work, which is based on combining two theoretical methods. First, we will 
adopt the quasi-chemical approximation (QCA) to the statistical mechanics of interstitial 
solutions [12,24-27]. In the QCA all interstitial sites are available for mixing, but the C 
atoms are regarded as exerting a repulsive force on each other, so that they enter adjacent 
interstitial positions less frequently as would be the case if their distribution were random. 
Thus the QCA will allow us to treat soft-blocking effects in austenite. The key parameters 
in this treatment are the energies of formation of the C-C and C-Va pairs, which will be 
accurately determined by analysing Oc data. Secondly, we will perform Monte Carlo (MC) 
simulations for various values of the Ncf^Pe ratio, using the pair formation energies from 
the quasi-chemical analysis. In this way, the average distribution of interstitials around a 
given C atom will be studied as function of composition. In particular, the average number 
of empty NNIS will be determined for various alloys, and compared with the results of 
previous studies, as well as with information extracted from Mossbauer experiments. 

In the analysis of the Mossbauer spectra of austenite various assumptions about the 
distribution of C in the octahedral interstitial sites have been proposed [28,29]. In particular, 
a model has been suggested for dilute solutions in which the 12 NNIS of a C atom are 
excluded [28]. In this case the C atoms occupy only the centre of the cubes of a structure 
with the formula FegCi-y, so that three possible environments for Fe atoms, associated to 
different hyperfine interaction may be distinguished, which are shown schematically in Fig.2, 
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viz., 

a. Fe atoms without nearest neighbor and next nearest neighbor C atoms, associated to 
the singlet Too- 

b. Fe atoms without nearest neighbour but with n next nearest neighbour C atoms 
(n = 1 — 4), ascribed to the singlet Fon- 

c. Fe atoms with one C atom nearest neighbour but without next nearest C neighbours, 
related to doublet Fiq. 

Alternatively, a model [29] has been proposed in which all the octahedral sites of the fee 
structure are available for occupation and no assumptions are made on the distribution of 
the C atoms in the second interstitial shell. The various Fe-C environments involved in this 
type of model are shown schematically in Fig. 3, viz., 

a. Fe atoms without nearest neighbours C atoms, associated to the singlet Tq. 

b. Fe atoms with one nearest neighbor C atom, or Fe atoms with two C atoms nearest 
neighbours at 90° from each other, related to the doublet Fi. 

c. Fe atoms with two C atoms placed at opposite nearest sites (180°), ascribed to the 
doublet 

The purpose of the present paper is to provide new experimental and theoretical infor- 
mation on the distribution of C atoms in the octahedral sites of the austenite phase. The 
work proceeds as follows. First, Mossbauer experiments are performed on a series of alloys 
with up to 12 at.%. 

II. EXPERIMENTAL 

A. Alloys, samples and heat-treatments 

Five Fe-C compacted graphite alloys were prepared in a medium frequency induction 
furnace using the sandwich technique in ladle to treat the liquid metal. The necessary 
amount of Si, Mg, Ce and Ca was added to obtain compacted graphite morphology. The 
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final alloy contained 3.40 wt.%C, 2.35 wt.%Si, 0.58 wt.%Mn, 0.04 wt.%Cu, 0.01 wt.%P, and 
0.02 wt.%S. 

Samples of 20 mm diameter and 3 mm thickness were taken from "Y-shape" blocks 
(ASTM A-395) cast in sand moulds. Samples SI, S2, S3, S4 and S5 (see below) were annealed 
at 1173K during 30 minutes, quenched in a salt bath and held at 623K during 1, 3, 4, 5 
and 10 minutes respectively. The samples for Mossbauer spectroscopy were prepared by 
conventional grinding techniques to reduce their thickness down to about 70mm, using 
diamond paste of 6, 1 and 0.1mm for final polishing. 

B. X-ray measurements 

X-ray measurements were performed in a Philips PW1710 diffractometer using the 
monochromatic Ka radiation of Cu, in Brag Brentano's geometry, with a step mode col- 
lection of 0.02, 10s by step, with 29 ranging from 39° to 98°. The X-ray patterns, presented 
in Fig.4, were analysed with the Rietveld method [30]. The actual C concentration in the 
samples was determined by combining the lattice-parameters (a) extracted from the diffrac- 
tion patterns with the known a versus composition relation for fee Fe-C alloys [28]. The 
resulting lattice-parameter and the inferred i/c values for the various alloys are listed in Table 
I. 

C. Mossbauer experiments 

Mossbauer spectra were taken in a transmission geometry using a blCoRh source of 
approximately 5mCi intensity and recorded in a standard 512 channels conventional constant 
acceleration spectrometer. In order to analyse in detail the austenite pattern, the spectra 
were taken in the velocity range between - 2 and -|- 2 mm/s. Velocity calibration was 
performed against a 12 mm thick a - Fe foil. All isomer shifts were referred to this standard 
at 298K. The spectra were fitted to Lorentzian line shapes using a non-linear least-squares 
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program with constraints. For the effective thickness of the samples analysed no Voigt 
line-shape correction was necessary [31]. 

The results of the Mossbauer experiments are shown in Fig. 5. The central subspectra 
were associated to austenite, whereas the external lines on the spectra to ferrite/martensite 
phases [28]. The hyperfine parameters and the relative fractions f;^ {l,m — number of C 
atoms in the first and second coordination shell, respectively) associated to the various Fe 
environments obtained using the models [28,29] referred to in Sect.I are listed in Table II. 
Concerning the second model [29], the contribution to the spectra of the doublet r2 associ- 
ated to Fe sites with two C atoms placed in opposite interstitial sites resulted undetectable 
by the present technique. 



The austenite interstitial solid solution is described as a lattice gas of carbon atoms 
and Nva vacancies, distributed in the N^-l- N^a = N^N^e octahedral interstitial sites of 
the fee structure (Fig. 1) associated to N^e iron atoms. The occupancy of the first and 
second interstitial shell was accounted for in three-dimensional MC simulations to calculate 
the number riij = C or Va) of C-C, C-Va and Va-Va pairs, the relative fractions f;^ 
associated to the different Fe environments, and the number Cjo of C atoms having i C 
atoms in the first interstitial coordination shell and none in the next interstitial shell. 

A Fortran 77 routine using the Monte Carlo method, an Ising-type Hamiltonian and 
periodic boundary conditions was developed. Metropolis method was used to define the 
probability of the C jumps. A randomly chosen C atom has the probability P to jump to 
an empty interstitial neighbouring site, also randomly chosen, viz.. 



III. THEORETICAL 



A. Monte Carlo simulations 




(1) 
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where Eti and £t/, are the initial and final total energies, respectively, calculated using the 
relation Et = nc-c^^- Here = 2 Sc-va - ^c-c is the energy of formation of a C-C 
pair of nearest neighbour C atoms relative to the individual C atoms (see below), nc-c is 
the number of C-C pairs, ec-c and ec-va are the interaction energies of the C-C and C-Va 
pairs, respectively, R is the gas constant and T is the temperature in Kelvin. If the atom 
movement decreases the total energy, the jump is allowed {P—1), but if the total energy 
increases, the jump is allowed with a probability P = exp[{eTi — eTf)/RT]. The Fe atoms 
remain still during the simulation, and their positions were only used to calculate the number 
of Uij pairs and the relative fractions f;^ associated to the different Fe environments. 

In order to study the convergence of the results, cells of 4^, 6^, 8^ and 10^ were used. 
For simulations using cell sizes of 6^ and higher the riij and iim fractions did not vary, hence 
cells of 864 Fe atoms and the corresponding number of C atoms were employed to decrease 
the calculation time. For all C concentrations, the equilibrium of the system was attained 
approximately at three MC steps, where a MC step is defined as Nc attempts of movement 
of a C atom. 

Finally, the occupation of the interstitial sites was characterised using the average number 
z of empty NNIS, which was calculated from the MC results as follows: 

EC,o(12-i) 

z-'^, (2) 

B. Quasichemical model calculations 

The energy of formation of a C-C pair that enters in the MC calculation was determined 
by analysing experimental Oc data in terms of the QCA to the statistical mechanics of the 
Fe-C solutions developed by Bhadeshia [27] . This formalism yields for the activity Cc 
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where Z (=12) is the number of NNIS, and AGc is the Gibbs energy of C in austenite 
relative to graphite. The value of the parameter A that minimise the Gibbs energy is: 

A = ^^{l-[l-4a|/e(l-|/c)^]} 



2(7 



with 



a — \ — exp 



A linear approximation of Eq.3, appropriate for describing the dilute solution range was 
fitted to carbon activity data measured at 1423K [22] . A least-squares fit of the equation: 

1 - 2/c 



i?T In 



Vc 



y^ZAe + AGc 



which is shown in Fig. 6 yielded AGc=445l25 cal/mol, and A£=149239 cal/mol. This one 
was adopted in the MC simulations. 

The number of pairs calculated for < yc < 1 using the quasichemical formalism 
(Table III) are plotted in Fig. 7 using lines. For comparison the Uij determined in the MC 
calculations using Ae = 1492 cal/mol are plotted using symbols. The inset there gives a 
comparison for the composition range corresponding to the experimental solubility of C in 
austenite, viz., t/c < 0.1. There is a very good agreement between the QCA and the MC 
predictions for n^j , which encourages a discussion of the MC results for the relative fractions 
iim of the various Fe environments, as functions of ye- 



IV. DISCUSSION 

A. Monte Carlo versus Mossbauer results 

The MC results for the lim as functions of yc are plotted in Fig. 8. According to the 
present simulations the main contributions to the Mossbauer spectra originate in the Fe 
environments without C atoms in the first interstitial shell (foo and fo™). Next in importance 
is the contribution of environments with one C atom in the first interstitial shell and none 
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in the second (fio) . Further, the simulations show that the contribution of Fe atoms having 
the first and second interstitial shells occupied (f„rn) is not negligible, which indicates that 
this kind of Fe environments should be accounted for in fitting Mossbauer spectra. Figure 
8 also demonstrates that the relative fraction associated to Fe atoms with more than one C 
atom in the first interstitial shell and without C atoms in the second shell (f„o) is negligible. 
This result contradicts the assumption of one of the models [29] developed to interpret the 
Mossbauer pattern of the austenite, which was reviewed in Sect. I. In fact, we find that the 
contribution of Fe environments with two C atoms either at 90*^ or 180° is negligible. 

In order to compare the MC results with the results of analysing the present Mossbauer 
spectra (Table II), the relative fractions foo and the sum fo = foo + fon were chosen as the key 
quantities in the FegCi-^ model [28] and the random model [29], respectively. The reason 
for this choice is that Fe environment with C atoms in the first and the second interstitial 
shells have not been considered by any of the models proposed to reproduce the Mossbauer 
spectra [28,29]. Moreover, we have also shown that the fraction of Fe environments with two 
C atoms in the first shell is negligible according to MC results. Hence, the only fractions 
that can be determined without ambiguity for these models are foo and fo, respectively. 

A comparison between the foo and fo versus composition values extracted from the 
Mossbauer spectra (open symbols) and the MC simulations (filled symbols) is presented 
in Fig. 9. According to Fig. 9 the "random model" underestimates significantly the con- 
tribution of the Fe environments without nearest neighbour C atoms. This suggests that in 
this composition range, some blocking effect of the NNIS in austenite should be accounted 
for. In agreement with this expectation, the MC results fit remarkably well with the foo 
values extracted from Mossbauer experiments when the FesCi-y model is adopted, i.e., the 
model including some blocking effect of the interstitial sites. 

Finally, MC simulations have been reported previously [32] which are based on extracting 
from Mossbauer measurements a weak C-C repulsion in the first coordination shell (wi = 
830 cal/mol) and a stronger one (w2 = 1730 cal/mol) in the second shell. Such calculations 
were interpreted [29,32] as indications that the FesCi-^ model is not adequate to represent 
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the austenite phase. Similar results (wi = 830 cal/mol) were arrived at in ref. [29] by 
assuming that pairs of C atoms can occupy the first interstitial shell either at 90° or at 
180°. The present MC results, based on energies extracted from thermodynamic data, do 
not support these results. 

B. Account of blocking effects 

In Fig. 10 the average number of empty NNIS z obtained from the MC simulations is 
plotted as function of yc (symbols). The z values corresponding to the composition of the 
present experimental alloys are plotted using empty symbols. The dashed line in this graphic 
refers to the z value corresponding to a random mixture, viz., z = Vl{\ — y^)- The empty 
symbols in Fig. 10 indicate that already in alloys with y^ — 0.05 the NNIS of the C atoms 
are, on the average, less occupied than in a random mixture. This fact is in qualitative 
agreement with the ideas behind the excluded-sites model, which motivates the following 
analysis of the blocking effects in models for Oc in austenite. 

It has recently been pointed out [23] that in a strict hard-blocking model, empty sites 
must be interpreted as blocked sites. This implies that the h parameter of the HBESM 
(Sect.I) should be considered as equal to z (Fig.lO). Two consequences of such interpre- 
tation will be discussed. The first consequence is that the MC results in Fig.lO cannot be 
represented using the HBESM unless the h parameter is allowed to vary with composition. 
In qualitative agreement with this. Gates et al. [21] interpreted their own z values from MC 
calculations as a composition dependent h parameter, decreasing with the increase in the C 
content. However, the excluded sites model does not explain the composition dependence 
of &, which has stimulated some attempts to improve the simple picture by invoking, e.g., 
overlapping of the sites excluded by different interstitial atoms [17,19,21,23]. The second 
consequence of the current [23] interpretation of the hard-blocking is that Fig. 10 yields h 
{—z) >10 in the composition range yc < 0.2. However, these values cannot be reconciled 
with those extracted from experimental Oc data. The latter are integral and non-integral h 
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values falling in the range 3 <b <5 [9,19,22,23]. 

In view of these facts we conclude that the strict form of the excluded-sites model (Sect. I) 
does not seem able to account for the present blocking effects by using the same b values 
which are known to reproduce the experimental Qc data. It is also evident that a more 
realistic account of such effects would require abandoning the one-parameter formula for Oc, 
i.e., what has been considered as the main advantage of the HBESM [23]. 

V. SUMMARY AND CONCLUDING REMARKS 

In the present study the energy parameters describing the interstitial solution of C in the 
fee phase of Fe have been obtained by analyzing experimental thermodynamic data in terms 
of the quasichemical approximation. These parameters have been used as input information 
in Monte Carlo simulations, and various key quantities have been obtained. In particular, the 
composition dependence of various fjj ratios, describing the relative weight of the various Fe 
configurations contributing to the Mossbauer spectra of Fe-C austenite, have been predicted 
and compared with those derived by modelling the Mossbauer spectra. In this way, two 
alternative models [28,29] for the contribution of the Fe environments have been tested. The 
present comparison between Mossbauer and theoretical results indicates that a description 
similar to the FcgC structure [28] should be preferred for the Fe-C austenite phase. Such a 
model [28] is usually associated to the blocking of some of the nearest neighbour interstitial 
sites by a C atom. However, the present results cannot be accounted for by the simplest 
hard-blocking excluded sites model, often used to provide a one-parameter formula for the 
activity of C in austenite. We believe that the picture of blocking effects in austenite 
emerging from the present study should be useful in further attempts to refine the current 
models for thermodynamics of interstitial solutions. 
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FIGURE CAPTIONS 

Figure 1: Interstitials sites of the austenite fee phase, a) ■ : C atom, b) • : Vacancy, c) 
dotted Une: C-Va pair, d) sohd Une: Va-Va pair and e) dashed hne: C-C pair. Fe atoms are 
placed in the corners of the cube and in the face centres. 

Figure 2: Fe environments in the FesCi-y model [28]. Filled and open circles correspond 
to C and Fe atoms, respectively. 

Figure 3: Fe environments in the random model [29] . Filled and open circles correspond 
to C and Fe atoms, respectively. 

Figure 4: X-Ray diffractograms of samples SI to S5. The bottom bar diagrams indicate 
from top to bottom: ferrite, austenite, martensite and C graphite. 

Figure 5: Mossbauer spectra recorded for samples SI to S5. 

Figure 6: Linear fit for activity data of ref.[22]. The values As — 149239 cal/mol and 
AGc = 445I25 cal/mol were determined. 

Figure 7: The number Uj^ (i, j = C or Va) of pairs C-C, C-Va and Va-Va. Squares, 
triangles and circles represent Va-Va, C-C and C-Va pairs respectively, obtained from Monte 
Carlo simulations using A£=1492 cal/mol. a) dash-dotted line: Va-Va, b) dashed line: C-C 
pairs and c) dotted line: C-Va pairs, calculated using the quasichemical model with the 
same Ae value. The inset gives a comparison for the composition range corresponding to 
the experimental solubility of C in austenite, viz., jc <0.1. 

Figure 8: The relative fractions f/m, (1, m = number of C atoms in the first and second 
coordination shell, respectively) associated to the various Fe environments obtained using 
Monte Carlo simulations, as functions of C content. 

Figure 9: The relative fractions foo (circles) and fo (diamonds) associated to the various 
Fe environments obtained from Mossbauer data (open symbols) using the models [28,29] 
referred to in Sect. I, compared with results from Monte Carlo simulations (filled symbols). 

Figure 10: The average number z of empty nearest neighbour interstitial sites of a C 
atom in austenite calculated as a function of composition using Monte Carlo simulations. 
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TABLE CAPTIONS 

Table 1: Lattice parameter (a) determined from the difFractograms using Rietveld [30] 
and corresponding to the samples SI, S2, S3, S4 and S5. The C content was determined 
using the empirical relation of ref.28. 

Table 11: Hyperfine parameters and relative fractions of the different Fe environments 
found in austenite using the models of refs. 28 and 29. 

Table III: Pair interaction energy and number of n^ pairs in the austenite phase obtained 
using the quasichemical model [27]. 
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Table III 
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Figure 1, K. Laneri, PRB 
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Figure 2, K. Laneri, PRB 




Figure 3, K. Laneri, PRB 



intensity (a.u) 




Figure 4, K. Laneri, PRB 




Figure 8, K. Laneri, PRB 




Figure 9, K. Laneri, PRB 




Figure 10, K. Laneri, PRB 



